Exact magnetization plateaus and phase transitions in spin-S" Heisenberg 
antiferromagnets in arbitrary dimensions 
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We generalize a class of Heisenberg antiferromagnets in one, two and three dimensions, which 
have been shown to exhibit magnetization plateaus for spin-^. In a certain parameter range of the 
general model, which is formally defined in D dimensions, we obtain the exact ground state(s) in the 
presence of an external magnetic field for arbitrary values of spin. In this range, the magnetization 
remains a constant as a function of the external field, except at some special values of the field 
where there is a jump from one plateau to the next. The plateaus are formed at certain specific 
fractions of the full magnetization which are determined by the spin and the lattice. Our general 
spin-5 result reproduces the known cases for spin-i in various lattices. Furthermore, we argue that 
outside the exact regime, the mechanism for the plateau formation is different. This results in first 
order phase transitions along some of the plateaus as the coupling constant is varied. We rigorously 
show the existence of such transitions for some particular cases. Finally, we numerically analyze a 
spin-1 model in one dimension using exact diagonalization to obtain its complete phase diagram. It 
agrees with our analytic results. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Kz, 75.45. +j 
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Many low-dimensional spin systems, especially in one dimension, respond discontinuously to the variations of an 
external magnetic fieldi. In such systems, magnetization curve forms plateaus as a function of the field. In most cases, 
plateaus are continuously connected to one another, whereas a more interesting scenario is when there are abrupt 
jumps between two plateaus. This usually happens at high fields in systems with localized magnon excitations^. 
There is a class of models for which magnetization plateaus of the second type are formed via a simple mechanism. 
The basic feature of these models is that they consist of units of spins coupled together such that the total spin of 
each unit is conserved^i^i^iSii. The first of its kind was introduced and studied at zero field by Gclfand^. The system 
consists of coupled dimers with the following Hamiltonian (see Fig. ^, 
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FIG. 1; The Gelfand ladder. The full lines represent the coupling J' and the dashed lines J. 

Note that the total spin on each dimer (5°, + S'jJ is conserved. A whole lot of eigenstates can be trivially 
constructed for the above Hamiltonian. Any state in which each dimer has a specific total spin and is maximally 
polarized in the 2-direction will be an eigenstate. Then, one expects the ground state to be a product of dimer-singlets 
for J' » J. A lower bound on J' such that the dimer-singlet is the ground state has been obtained by rewriting the 
Hamiltonian as sum of units consisting of three spinal, and is given by. 
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For S = ^, the ground state in fact remains to be the dimer-singlet for J' > 1.4 J. For J' < 1.4 J, all the dimers go 
into triplets and form the ground state of a spin-1 uniform chain. 



The lower-bound for J' given in ^ suggests that, the higher the spin, the stronger the dimer bond needs to be for 
the dimer-singlet to be the ground state. But it turns out that the dimer- singlet is quite robust even for very large 
spins. In this paper, we show that the lower bound on J' can be pushed down to 2 J for arbitrary values of spin, i.e., 
the bound given by |(2Jl for S' = ^ is in fact true for all values of S. 

Honecker, Mila and Troyes^ (from hereon referred to as HMT) have numerically studied the 5 = ^ system in the 
presence of magnetic field. For J' > 2 J, the system starts with the dimer-singlet ground state at zero field and the 
magnetization per site (M) is zero. The system remains in the dimer-singlet phase till a transition field Bc^ = J', at 
which one set of alternate dimers form triplets polarized in the direction of the magnetic field (M =2)- For further 
increase in the field strength, this state remains to be the ground state till a second transition field Bc2 = 2 J + J', 
at which the rest of the dimers also become triplets and the system is fully polarized {M = 1). This behaviour has 
also been inferred on the basis of a strong coupling analysis^. For J' < 2 J, gapless phases come into play and the 
magnetization curve starts developing tails between plateaus before the latter disappear altogether. HMT also studied 
a three-leg ladder and obtained abrupt jumps between plateaus in some range of the parameter. Similar results for 
S = ■^ have been obtained for a modified Shastry-Sutherland model in two dimensions^ and a three-dimensional 
version of Gelfand ladder—. 

In this paper, we define a general model in arbitrary dimensions which includes all the models mentioned above. In 
a certain regime of the parameter space, we obtain the exact ground state(s) in presence of magnetic field for arbitrary 
values of spin. The magnetization curve forms plateaus and they jump abruptly from one to another at certain special 
values of the field. The magnetization on the plateaus are some rational fractions of the full magnetization, which 
depend on the underlying spin and the lattice. 

Outside the exact regime, we argue that the plateaus will survive in the immediate neighbourhood, quite possibly 
developing tails between successive ones. But the mechanism for plateau formation in this case is different from the 
exact regime and this results in first order phase transitions along some of the plateaus. We rigorously show the 
existence of such transitions for some special cases. Finally, we numerically analyse the ladder defined in J^l for 5 = 1 
using exact diagonalization. We obtain the complete phase diagram which confirms our analytic results. 

II. THE GENERAL MODEL AND THE EXACT GROUND-STATES 

HMT have defined a general n-leg ladder in one dimension. We generalize further and formally define Hamiltonians 
on an arbitrary bipartite lattice in D dimensions, where a bunch of n spins live on each lattice-site and interact 
among each other as well as with the spins living on neighbouring sites. The reason for writing down the most general 
Hamiltonian is that, most of our analysis is independent of details like dimension, lattice and spin. 

Let Sp.x (/^ = li 2, ■ • ■ , n) be the spins living at the site x. Define, 

n 

Tx = ^ S,,,x. (3) 

Then the Hamiltonian is defined as, 

^=yETx + ^ETx•Ty-i3^^x^ (4) 

X (x.y) X 

where ( , ) denote nearest neighbours. We take both the couplings to be positive. The models mentioned previously 
are particular cases of the above Hamiltonian. For the Gelfand ladder and the three-leg ladder in HMT, the underlying 
lattice is a simple chain and n = 2 and 3 respectively. Both for the modified Shastry-Sutherland model^ and the 3-d 
model of Sutherland^ n = 2. For the former the underlying lattice is square and for the latter it is one- fourth depleted 
cubic. 

The variables T^^s are also SU{2) generators, and they can take the representations having the Casimir 
Tx = j{j + 1); such that j = 0,1, ■ ■ ■ , nS, when nS is an integer; and j = i, |, ■ • ■ , nS, when nS is a half-odd integer. 
Thus the bunch of spins at each site kinematically act respectively like a bosonic or fermionic rotor, except that the 
values j can take are truncated at nS. The Hamiltonian in Eq. Q, feels the value of the underlying spin only through 
this truncation. This suggests that one can treat the problem for general values of spin. 

Now we will use a method of divide- and- conquer to solve for the ground state. If one can write H as a, sum over 
smaller units hi, such that there exists a state |^o) which is simultaneously a ground state of each and every hi, then 
^/jq) will be a ground state of H. To this end, we rewrite H as a. sum over all the bonds of the underlying lattice. Let 






FIG. 2: Writing the ladder Hamiltonian as a sum over units of two diniers. 



i label the bonds and let z be the co-ordination number of the lattice. Then 
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where Ti_i and T2,i are the two degrees of freedom living on the bond i. Note that the same spin will be shared by 
many bonds and hence labeled in more ways than one. For the Gelfand ladder, the above break-up is schematically 
represented in Fig [3 

Now we will concentrate on a particular hi, and find its ground state and energy as a function of B. For notational 
convenience we shall drop the index i from hereon. After a rearrangement of terms, the bond Hamiltonian can then 
be rewritten as. 



/i = ^ (J' - zJ) (T? 
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The four terms in h mutually commute and therefore the spectrum is trivially solved. Let ji(ji -I- 1), J2(j2 + 1), 
j{j + 1) and TO be the eigenvalues of Tj, T^, (Ti + Ta)^ and {Tf + r|) respectively. Here |(ji - J2)\ <j< (ji -|- J2) 
and —j < m < j. Then the spectrum for h is given by, 

E{ji,j2,j,m) = — (J' - zJ) (ji(ji + 1) + J2(J2 + 1)) + 2 Jij + 1) - - m. 



(7) 

Since the magnetic field couples to a conserved quantity, the problem of finding the ground state at a given field 
strength reduces to finding the lowest energy state in each j-sector at zero field, i.e., we need to minimize the field- 
independent part of E with respect to ji and 72- For non-zero field, the degeneracy arising from the SU{2) symmetry 
will be lifted and the ground state will have to = j. Next we consider the case J' > zJ, for which the coefficient of 
the first term in E is positive. 

When nS is an integer, ji, J2 = 0, 1, • ■ • , nS and for half-odd integer values, ji,J2 = 51 f''' ' ' '^'^- -'^^ both the cases, 
j = 0, 1, • • • , 2nS. Anticipating the solutions, we divide all possible cases into two: 

1. Even j, integer nS; 

Odd j, half-odd integer nS. 

2. Odd j, integer nS; 

Even j, half-odd integer nS. 

For Case 1, E is minimized by the choice. 



Ji = J2 



And for Case 2, the minimum of E is when. 



Ji = 



(.1 - 1) 



J2 



(j + 1) 



(8) 



(9) 



There is one exception to the above general rule, though. When nS is half-odd integer and j ~ 0, Eqs. © give 
ji = — i and J2 = h- This is not an admissible solution and the correct solution is, ji = J2 = -^^ We will have more 
to say about this exceptional case later. 

The ground state energies in the presence of field corresponding to the solutions (jH)) and are respectively given 
by, 
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Let \l,m)a denote the state with eigenvalues l{l + 1) and m for the operators T^ and T^ respectively. Then the 
corresponding ground states are, 

\4'\j)) = IjA J72)i®b72, j72)2, (12) 

\4'\j)) = I(J - l)/2, (.7 - l)/2)i ® |(j + l)/2, (j + l)/2)2. (13) 

(12) 

Since the field-independent part oi Eq are monotonically increasing functions of j, as B is increased, j will go up 
sequentially in steps. The value of the field at which the total spin goes up to (j + 1) from j can be obtained from 
Eqs. lfTn|l and ((TT|l . For Case 1, the transition field is, 

B^ = ^-^^ 3 + J'. (14) 

For Case 2, it is given by, 

Bif=^-^^^^{j + l). (15) 

This completes the solution of the bond-Hamiltonian {h) for J' > zJ . The important point to note is that the ground 
states ()12|l and (|13|l are product states of the two sites connected by the bond. In the context of the full lattice, this 
means that we can put together these states living on the bonds and construct a state which will be a simultaneous 
ground state of hi for every i. Consequently, it will be a ground state of the full Hamiltonian, H . 

For Case 1, the ground state consists of all dimers having total spin j/2 and maximal polarization in the z-direction, 
i.e., T^ = (j/2)(j/2 + 1) and T^ = j/2. (Note that j denotes the total spin on the bonds, which are not good quantum 
numbers. Here we merely use j to label the ground states). Since the ground state is unique for individual bond 
Hamiltonians, this state has to be the unique ground state of the full Hamiltonian. The state can be written as, 

l*o'^(j))=ni^V2,j72)x (16) 

X 

For Case 2, the ground state can be written as, 

l*g(j)) = n 10' " l)/2' 0^ - l)/2)x n I(J' + l)/2' 0' + l)/2)y. (17) 

xGA yG-B 

where A and B denote the two sub-lattices. The above ground state breaks the translation symmetry and is therefore 
doubly degenerate, the other ground state being the one in which the two sub-lattices are interchanged. Eqs. I|16() 
and H17|) give the ground states for all values of j ranging from to 2nS. The fractions of the total magnetization at 
which the plateaus form are, 

and the transition fields are given by Eqs. H14() and H15|l . At each jump of the magnetization, dimers belonging to 
one of the sub-lattices alternately get excited to a total spin one unit higher. These plateau transitions are similar to 
the macroscopic magnetization jumps discussed by Schulenburg et. al^. The difference being that, in our case the 
localization of magnons arises from conservation laws. 

The above solution for the ground state hinges on one crucial aspect - the ground state of the bond Hamiltonian 
(h) is a product state on the two sites. For half-odd integer nS and j = 0, we have seen that the ground state of 
h has ji — J2 — \- This is not a product state and thus cannot necessarily be satisfied for all bonds. When nS is 
half-odd integer, the minimum spin on any site is ^ and therefore at zero field the system will be in the ground state 
of spin- 2 system on the underlying lattice. This is gapless, even in one-dimension^, and therefore the magnetization 
will vary continuously with the field. For j = 1, the solution again has ji = J2 = \- This is the fully polarized state of 
the spin- 2 system. The field at which it becomes fully polarized can be determined by the gap of one-magnon states 
(i.e., all spins except one pointing up). For a hyper-cubic lattice in D dimensions, this is calculated to be, 

^y = ASD (19) 

Our general expressions for the magnetization fractions (Eq. (|18|l ') and transition fields (Eqs. (|14|l . (|15|l and (|19|l ') 
check correctly with the results in HMT for spin-i, two and three-leg ladders. 



A comment is in order here regarding the ground states at plateau transitions. Suppose the field is tuned to Be , 
i.e. the transition field at which j goes up to j + 1. Also let j be even (a similar argument holds for odd j as well). 
Then the state of each bond can be any one of the following three states. 

\MJ))i = li/2, j/2) ® \j/2, J 12), (20) 

\M3))2 = U/2, j72)® 1072 + 1), 072 + 1)), (21) 

\Mj)h - 1072 + 1), (j72 + l))®b72, .7/2). (22) 

This gives rise to a degree of degeneracy which grows exponentially with the system size. In HMT, it is mentioned 
that the ground state of the spin-^ ladder is either all triplets (not necessarily polarized) or a mix of singlets and 
fully polarized triplets. They make an exception for those special values of the field at which there is a jump in 
magnetization. The above analysis shows that even at transition fields the ground states belong to the afore-mentioned 
set. 

III. BEYOND THE EXACT REGIME 

The divide- and- conquer method we employed in the previous section to find the exact ground state worked because 
the ground state of the bond Hamiltonian [h) is a product state. This is no longer true for J' < zJ. Then the ground 
state of h for a particular j will have, 

Tl^nS{nS+l), Tl=nS{nS+l), (Ti + Ta)^ - jO + 1), (^f + r|) = j, (23) 

This state is not a product state unless j = nS. But this does not necessarily mean that the product states cease to 
be the ground state of H for all J' < zJ . Now we will argue that the product states, and hence the magnetization 
plateaus, do indeed survive for a while as J' goes below zJ . Significantly, the solutions are not same as that for 
J' > zJ as given by Eqs. (|16|l and H17(l . The new solutions also form plateaus at the same fractions and it results in 
a first order phase transition on each plateau (with certain exceptions) at J' — zJ . 

For the proof of the existence of product state solutions for J' < zJ, we do the following. First we will show 
that, at J' = zJ , the only solutions are a finite number of product states, except at the transition fields. Then, by 
assuming these states to be gapped, it will follow that there exists a positive non-zero A such that the ground state 
at (z — A) J < J' < z J will be among the ground states at J' = zJ . 

At J' = zJ , the bond Hamiltonian in Eq. ^ becomes, 

/i = ^(Ti+T2)2-|(Ti^+T|). (24) 

The above Hamiltonian sees only the total spin of the bond, but is independent of the spins of the sites. Then the 
ground states of a given j-sector are highly degenerate - there is a freedom to choose ji and j2 as long as they combine 
to form total spin j. Among these states, there is a subset given by, 

|^o)p = |^7)i® 10 -0,0-0)2, / = 0,l,---,j, (25) 

which constitutes all the only product states. Next we will show that the ground states of the full system at J' = zJ 
can be written as product states on the sites. Let us consider two neighbouring bonds - a connecting sites 1 and 2, 
and h connecting sites 2 and 3. Also let ha and /ib be the corresponding bond Hamiltonians. Then, at J' = zJ any 
ground state, of {ha + /»&) has to satisfy the following conditions. 

(Ti+T2)'|0o) = jO + 1)I0o), (26) 

(27) 

l)l0o), (28) 

(29) 

where j is determined by the field, B. In the sub-space of fixed eigenvalues for T^, Tj and Tg, Eqs. H2()|) and 127|) 
imply that, 

|0o) = |0o)l2 ® |0o)3, (30) 



(TJ + T^)|0o) 


= i\<t>o), 


(T2 + T3)'|0o) 


= i(j + 


(T^ + T^)|0o) 


= j|0o), 



where, \(j)o)i2 lives in the combined Hilbert space of the sites 1 and 2, and |</>o)3 hves in the Hilbert space of the site 
3. This follows from the fact that for fixed values of ji and J2, there is a unique state in the Hilbert space of sites 1 
and 2 which satisfies Eqs. l(^ and (|?7jl . Similarly Eqs. ((SHJ and (|^ imply that, 

|0O> = \Ml ® l'/'0)23, (31) 

where, \(f>o)i lives in the Hilbert space of the site 1 and \4>o)23 lives in the combined Hilbert space of the sites 2 and 
3. Finally, Eqs. ^ and (EH imply that, 

l^o) = \h)i ® \M2 «> IMi, (32) 

i.e., |(/>o) is a product state of the three sites. Since every site is part of at least two bonds, this necessarily implies 
that the ground states of iJ at J' = zJ are product states. These ground states can be written as, 

\^oiJ'^zJ,j))^l[\l,l)^®Y[\j-l,j-l)y, / = 0,l,---,j, (33) 

xeA yes 

What we have shown now is that, though the bond Hamiltonian (/i) admits solutions which are not product states at 
J' = zJ, such solutions are not admissible for the full Hamiltonian (H). 

At transition fields, there is a greater degree of degeneracy. Then the constraint on the ground state is more relaxed 
- the total spin on a bond can be j or j + 1. This gives rise to a degeneracy which is local and hence grows exponentially 
with the system size. Away from the transitions fields, the ground states are those given by Eq. H33|l and are not 
related to each other by any local transformation. It is then fair to assume that there is a gap to excitations. 

Now let us see why gap at J' = J implies the existence of product ground state for {z — A)J < J' < 2 J, for some 
positive A. The three terms in H in Eq. |0J mutually commute. This means that one can choose the eigenstates 
to be independent of the three parameters J', J and B. The energy spectrum will depend linearly on all the three 
parameters and all transitions of the ground state are via level crossings. (The transitions will be first order unless 
there is a continuum of level crossings - a possibility in the thermodynamic limit). A gap to the ground states at 
J' = zJ would then imply that there exists some non-zero positive number A, such that for (z — A) J < J' < zJ the 
ground state(s) will be among those at J' = zJ. Among the states in Eq. 133|l . the one with the lowest energy for 
J' < zJ is given by. 



ieA keB 

l*o<(j)) = n |2^'2^) ® n K^' - 25), (j - 25)), for j > 2S. (35) 

ieA keB 

These states also have the same fractions for magnetization as before. But, now the system prefers to put all the spin 
into sites in one sub-lattice, if that is possible. This is in contrast to the situation for J' > zJ where it is energetically 
favourable to distribute the spin equally among the two sub-lattices. At j = 25*, one sub-lattice will be saturated and 
then the other sub-lattice start getting excited to higher and higher spins. 

In the above states, when j > 25, sites in one sub-lattice have spin 25* and those in the other sub-lattice have 
spin (j — 2S). On the plateau, each spin is maximally polarized in the z-direction. Below a critical field, gapless 
phase of the same set of spins will have lower energy. The critical field, in units of J, will be the gap in the spectrum 
of onc-magnons (states with all spins but one maximally polarized, the odd one having z-component one less than 
maximum). For a hyper-cubic lattice in £)-dimensions consisting of spins j'l in one sub-lattice and j2 in the other, 
the critical field is given by, 

^=2Din+j2). (36) 

HMT have numerically found such lines of transition for S — ^, two and three legged ladders and the formula above 
is consistent with their results. For the two-leg ladder, the transition occurs when the fully polarized effective spin-1 
chain becomes gapless. Then ji = j 2 = 1 , and Bcrit = 4 J. For the three- leg ladder there are three cases - a) jx = i, 
32 — \, b) ji = ^, J2 = I and c) ji = |, ^2 = §■ Then the transitions occur respectively at Bcrit = 2 J, 4 J and 6 J. 

To summarize our results thus far, we first solved for the ground states of the general spin-S' model in arbitrary 
dimensions for J' > zJ. The product state solutions we obtain exhibit magnetization plateaus. Then we argued that 
there exist product state solutions and hence magnetization plateaus even as J' goes below 2; J, but the mechanism 
for plateau formation is different from the case J' > zJ. Thus, on a plateau, there are first order phase transitions at 
J' = zJ. 
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FIG. 3: Writing the Hamiltonian as a sum over units of three dimers. Bold fuU hues represent bonds of strength J' , thin full 
lines J'/2 and dashed lines J. 



HMT have obtained a similar first order phase transition on the plateau for the three-leg, S = ^ model. 



In their 
case the transition occurs on M = ^ plateau from a fully polarized S ~ ^ chain to an alternating chain consisting of 
spins 5 and |. The latter is not a simple product state. 

The plateau solutions for J' > 2J, given by Eqs. lfTC|) and ((T7|l . and that for J' < 2J, given by Eqs. l|M|) and (jSHl, 
coincide for j = 0, 1, {2nS — 1) and 2nS. For S ~ ^ ladder {n = 2), all the plateaus belong to the above set and 
therefore there is no first order transition at J' = 2 J along any of the plateaus. 

From now we will concentrate on the two-leg ladder. First we will show for some higher spin cases that first order 
transitions along the plateaus at J' = 2J does indeed exist. 



IV. PHASE TRANSITIONS FOR THE LADDER AT J' = 2J 



We have seen that the ground states of the bond Hamiltonian are not product states for J' < 2 J. Therefore, to 
look for product state solutions in this range, we try a different divide- and- conquer scheme involving bigger units. To 
write Hi as sum over units containing three dimers then suggests itself ( see Fig. OJ, 



H, 



^j:^. 



(37) 
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JT2. • (T2,_i + T2,+i) - I (r|,_i + 2 T^, + T^VO . (38) 



As before, we drop the index i and write the 3-dinier Hamiltonian as, 



^-t(t? 



2T? 



JT2 • (Ti + T3 



|(Tf + 2T|+T|). 



(39) 



T^, T|, T§, (Ti -I- Ta)^ and (Tf -I- T| -I- T|) are mutually commuting conserved quantities and this essentially reduces 
the problem to that of two (variable) interacting spins coupled to a magnetic field of unequal strengths at the two 
sites. Unfortunately, the spectrum cannot be found analytically for arbitrary values of spin. 

We have numerically diagonalized h for spins up to 3. We have found that the product state solutions given in 
Eqs. (|34|l and H35|) do indeed exist for some J' < 2J. (This diagonalization can be easily performed for higher spins 
also, but we chose to restrict it to reasonable values). Note that the existence of the product state solutions for h 
necessarily implies their existence for the full ladder. This in turn confirms our prediction of first order transitions on 
a plateau at J' = 2 J. 

Next we study the particular case S = 1. We numerically analyse it to obtain the phase diagram in the complete 
parameter space and thus confirm our analytic results. 



V. SPIN-1 LADDER 



The spin-i ladder was numerically studied in great detail by HMT to obtain the complete phase diagram. Our 
results for J' > 2 J and the arguments for the existence of the plateaus beyond the exact regime are consistent with 
their phase diagram. As we mentioned in the previous section, to obtain the new phase predicted by us, one needs to 
look at spins higher than ^. In this section we present the results of an exact diagonalization analysis of the spin-1 
ladder. 

Since the individual dimer spins are conserved, we need to analyse only effective spin chains consisting of spins of 
varying values at each site. We have employed the Lanczos algorithm to exactly diagonalize the Hamiltonian and 
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FIG. 4: Phase diagram for the spin-1 ladder. Thick and thin lines denote first and second order transitions respectively. 




FIG. 5: Magnetization curves for the spin-1 ladder. The fractional values of magnetization are denoted on the plateaus. 



obtain the ground state energies. We have worked with 10 diniers (20 spins) and used periodic boundary conditions. 
For the finite system that we studied, we have checked that all the relevant states are such that the spins of the dimers 
belonging to a given sub-lattice are the same. There arc exceptions to this at plateau transitions, arising from the 
local degeneracy which has been discussed in Section |n] In the phase diagram (see Fig. 0)), [ji — J2]p and [j'l — J2]g 
respectively denote the fully polarized and the gapless phases of an alternating chain consisting of spins ji and J2 ■ In 
Fig. [SI the magnetization is plotted as a function of B/J for various values of J' /J. 

For J' > 2 J the analytic results from Section^ are very accurately reproduced. Since all the relevant states are 
product states, there are no finite-size corrections and agreement with the analytic results is exact up to numerical 
inaccuracies. There are five plateaus at the fractions 0, 4, ^, j and 1. The successive transitions occur at the 



following values of the field, 
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(40) 



:^ = 2 + 2-, :^=4 + 2-. (41) 

J J' J J ^ ^ 

On the i\/ = i plateau, [1 — l]p phase gives way to [2 — 0]p at J' = 2 J, as we expect from our analysis in Section ITTll 
From [2 — l]p and [2 — 2]p phases there are second order transitions to the corresponding gapless phases [2 — l]g 
and [2 — 2]g respectively. The lines of these transitions, obtained from Eq. (|36|l . are given by Bcrit = 6J and 8J 
respectively. These lines meet the lines separating plateaus, given in Eq. (|41() at J' — 2 J. 

For J' ^ 2 J, all the dimers have spin 2. Then the ground state at zero field is the Haldane phase of the uniform 
spin-2 chainiiii^. When the field is turned on, the magnetization remains zero till a critical field at which the system 
becomes gaplessi^. Then with further increase in field, the magnetization continuously increases and finally reaches 
saturation ([2 — 2]p). In the plot for magnetization curves in FiglSj there are discernible steps even in the gapless 
phases [2 — l]g and [2 — 2]g. These arc spurious jumps arising from finite-size effects. To obtain more accurate results, 
it will be necessary to employ techniques like DMRG. 

In the intermediate range, i.e. around J' ~ 1.6 J, just like in the S = ^ case in HMT, with changing field, the 
system goes in and comes out of various gapped and gapless phases. At zero field, [0 — 0] phase goes into spin-2 
Haldane phase at J' « 1.5975J. 

For higher spins, there will be (45 H- 1) plateaus for J' > 2 J as discussed in Section^ First order transitions of the 
type, [1 — l]p to [2 — 2]p will then happen at all plateaus except the first two and the last two. There will also be second 
order transitions of the type [2 - l]p to [2 - l]g, between [2S - (j-2S')]p and [25' - {j-2S)]g for (25-t-l) < j < 45. 
Thus we expect most features of the phase diagram in Fig. 01 to be generic for all spins. 

VI. SUMMARY AND CONCLUSION 

We have generalized a class of existing spin-2 models in various dimensions, which exhibit magnetization plateaus, 
to arbitrary values of spin. By using simple techniques we are able to analytically solve the spin-5 model in a 
substantial regime of the parameter space. Our approach also throws light upon the phases beyond the exact regime. 
In particular, we are able to predict that there will be a new type of plateau phase. Such phases cannot exist for 
the spin-i ladders studied in HMT. We have tested our results numerically for the particular case of 5 = 1 two-leg 
ladder. 

Though our general model is defined on arbitrary bipartite lattices, each site has a sub-structure which includes 
n spins. Therefore, constructing physically feasible models, where bonds of equal strength have equal length, is a 
non-trivial task. As far as we know, the one-dimensional linked tetrahedra introduced by Gelfand^ and Sutherland's 
generalization of the same to 3-d are the only possibilities. 
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